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Using first-principles calculations, we explore the electronic and magnetic properties of graphene 
nanomesh (GNM), a regular network of large vacancies, produced either by lithography or nanoim- 
print. When removing an equal number of A and B sites of the graphene bipartite lattice, the 
nanomesh made mostly of zigzag (armchair) type edges exhibit antiferromagnetic (spin unpolarized) 
states. In contrast, in situation of sublattice symmetry breaking, stable ferri(o) magnetic states are 
obtained. For hydrogen-passivated nanomesh, the formation energy is dramatically decreased, and 
ground state is found to strongly depend on the vacancies shape and size. For triangular shaped 
holes, the obtained net magnetic moments increase with the number difference of removed A and 
B sites in agreement with Lieb's theorem for even A+B. For odd A+B triangular meshes and all 
cases of non-triangular nanomeshes including the one with even A+B, Lieb's theorem does not hold 
anymore which can be partially attributed to introduction of armchair edges. In addition, large 
triangular shaped GNM could be as robust as non-triangular GNMs, providing possible solution 
to overcome one of crucial challenges for the sp-magnetism. Finally, significant exchange splitting 
values as large as 0.5 eV can be obtained for highly asymmetric structures evidencing the po- 
tential of GNM for room temperature carbon based spintronics. These results demonstrate that a 
turn from 0-dimensional graphene nanoflakes throughout 1-dimensional graphene nanoribbons with 
zigzag edges to GNM breaks localization of unpaired electrons and provides deviation from the 
rules based on Lieb's theorem. Such delocalization of the electrons leads the switch of the ground 
state of system from antiferromagnetic narrow gap insulator discussed for graphene nanoribons to 
ferromagnetic or nonmagnetic metal. 

PACS numbers: 81.05.ue, 75.70.-i, 75.75.-c 



I. INTRODUCTION 



Two-dimensional graphene has emerged as a natural 
candidate for developing "beyond CMOS" nanoelectron- 
ic^"iS. In addition to the reported huge charge mobilities, 
the weak intrinsic spin-orbit coupling in carbon-based 
sp^ structuresii^ could potentially allow for very large 
(micron long) spin diffusion lengths. These features, to- 
gether with the other " semi-conductor like" properties of 
graphene, make graphene-based spintronic devices highly 
promising^iiS and have triggered a quest for controlling 
spin injection in grapheneii"— . Many routes have been 
attempted to induce magnetism by proximity effect or 
inject spins from magnetic electrodes^^. Another, more 
intrinsic, possibility is shaping the geometry of graphene 
by designing graphene nanoribbons with zigzag edges. 
This has been found to induce localized edge magnetic 
states which can serve as a conceptually new building 
block for spintronic^—. 



The existence of intrinsic magnetism driven by atomic- 
scale defects (such as vacancies, chemisorbed species, 
grain boundaries,...) has been also suggested theoreti- 
callyi^i^ — , but remains fiercely debated on the experi- 
mental side^. It is indeed particularly difficult to achieve 
a precise experimental characterization of those defects, 
whereas the control of their density, positioning, or chem- 
ical reactivity seems an insurmountable challenge, jeop- 
ardizing a further use of magnetic properties in real de- 
vices. Additionally, the absence of a true energy gap in 
two-dimensional graphene limits the elaboration of ac- 
tive graphene-based devices and circuits with standard 
semiconductor technologies. 

Other route to make graphene magnetic is either 
chemisorption of odd number of adatoms or functional 
group^i^, or using magnetism on zigzag edge^*^!^. 
In the first case the stability of magnetic configurations 
at room temperature can be easily destroyed by the mi- 
gration of adatoms with turning the system into nonmag- 
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netic configuratior>2S. In contrast to the adatom based 
magnetism edge^ and vacanc}*^ magnetism in graphite 
is stable at room temperature. But herewith locahzation 
of the magnetic moments on the edges provides formation 
of the AFM exchange interactions between two edges^^. 

The case of graphene nanoribbons obey Lieb's theorem 
because the locahzed electrons on one edge belong to sub- 
lattice A and localized electrons from other edge to sub- 
lattice B. Magnetism on the edges of graphene nanoflakes 
is also described by this theorem^. Electron localization 
plays an important role in the different many-body effects 
on graphene edges^^ and bulk graphene^^. Understand- 
ing the nature of the electron localization and delocal- 
ization in graphene and related systems is necessary not 
only for control and manipulation of magnetism in stud- 
ied compounds but also for the development of knowledge 
about systems with strongly correlated electrons. 

Graphene nanomesh (further GNM) is the intermedi- 
ate compound between graphene nanoribbons with local- 
ized electrons on zigzag edges and perfect bulk graphene 
with delocalized electrons. The fabrication of GNM, us- 
ing block copolymer lithography and offering versatil- 
ity in varying periodicities and neck widths down to 5 
nm^, could circumvent the hurdles. Indeed, such tech- 
nique allows a scalable engineering of superlattices of 
large graphene vacancies whose density, shape and distri- 
bution can be controlled down to the nanoscale. Addi- 
tionally, GNM-based field-effect transistors were shown 
to withstand current densities two orders of magnitudes 
larger than individual graphene nanoribbon devices, with 
comparable on/off ratio and easily tunable by varying 
the neck width. Moreover, in accordance to the well- 
established Lieb's theorem^^, the possibility to control 
inner edge structures'^ of nanomesh vacancies could en- 
able a true control of intrinsic magnetic properties. 

In this paper, we present first-principles calcula- 
tions of electronic and magnetic properties of graphene 
nanomesh, and found that by varying the shape, differ- 
ent types of intrinsic ferrimagnetic'^ states can be ob- 
tained with clear identification of optimal conditions. 
Systematic studies of non-passivated and hydrogen- 
passivated GNM are achieved by varying the difference 
{Aab = 1^ ~ ^1) between missing A and B sites of the 
underlying bipartite lattice and analyzing different hole 
geometries. For non-passivated GNM with Aab = 0, 
stable non-magnetic states are found for armchair edge 
termination, while zigzag edge terminations result in 
antiferromagnetic ground states. These localized mag- 
netic moments however vanish when all edge C- atoms are 
hydrogen-passivated. In sharp contrast, when Aab 7^ 0, 
stable ferri(o) magnetic states are induced with net mo- 
ment up to 4 /iB (per 6x6 unit cell) originating from 
dangling bonds of edge atoms. 

Furthermore, for hydrogen-passivated GNM, the for- 
mation energy is dramatically decreased, and ground 
state is found to strongly depend on the vacancies shape 
and size. Our calculations reveal the existence of three 
magnetic regimes which depend on Aab- (i) highly mag- 



netic GNMs obeying Lieb's theorem corresponding to tri- 
angular shaped holes with even A+B; (ii) GNMs with 
quenched magnetic state due to complete chemical bond 
reconstruction with Aab = 1 and trivially nonmagnetic 
state with Aab = 0; and (iii) GNMs following interme- 
diate regime between magnetic and quenched magnetic 
states, i.e. triangular GNMs with odd A+B and more 
complicated structures including both even (e.g. sector 
shaped GNM) and odd (e.g. pentagon shaped GNM) 
A+B. We show that large triangular GNMs could be 
as robust as non-triangular GNMs providing possible so- 
lution to overcome one of crucial challenges for the sp- 
magnetism. Moreover, significant exchange splitting val- 
ues as large as ~ 0.5 eV can be obtained for highly asym- 
metric structures evidencing the potential of graphene 
nanomesh for room temperature carbon based spintron- 
ics. 



II. METHODS 

First-principles calculations were performed using Vi- 
enna ab initio simulation package (VASP)^ based 
on density functional theory (DFT) with generalized 
gradient approximation (GGA) for exchange correla- 
tion potential. We have used projected augmented 
wave method (PAW) — with the Perdew-Becke-Erzenhof 
(PBE) parametrization"^^ potentials to describe the core 
electrons of carbon. Periodic 6 by 6 unit cells were used 
to simulate non-passivated GNM structures as shown in 
Fig. [11 whereas periodic 8 by 8 unit cells were used to sim- 
ulate H-passivated GNM structures. The kinetic cutoff 
energies for the plane wave basis set used to expand the 
Kohn-Sham orbitals were 520 eV for the self-consistent 
energy calculations. Methfessel-Paxton method^^ is used 
with a broading width of 0.2 eV for the partial occupi- 
cancy smearing calculations. A 9x9x1 k-point mesh was 
sufficient to ensure good convergence in the total energy 
differences. The structural relaxations were performed 
ensuring that the Hellmann-Feynman forces acting on 
ions were less than 10~^ eV/A. 



III. MODEL OF GRAPHENE NANOMESH 

A GNM can be formed by either removing atoms cen- 
tered of a six ring structure [Fig. [Ua,b,c)] or centered of 
one carbon atom [Fig. [Hd)], which either lead to form 
GNM with balanced or unbalanced number of removed 
A and B sites. For the sake of clarity, we label those 
structures according to their shapes and put the num- 
ber of removed A and B site atoms as a subscript. For 
instance, the structure of Fig. [TJa) named 63:3, corre- 
sponds to a Circle hole shape GNM with 3:3 denoting 
3 A and 3 B atoms removed from perfect graphene. The 
superscript H is used for hydrogen-passivated GNM. 
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FIG. 1: (a)-(c): Schematics of the calculated crystalline struc- 
tures for balanced non-passivated circular shaped €3:3, C6:6 
and Ci2:i2 GNM structures, respectively; (d) the same for un- 
balanced non-passivated triangular shaped Te:? GNM struc- 
ture. Edge carbon atoms are in blue and red color to repre- 
sent A and B sites, respectively. For convenience, positions 
of removed atoms are indicated in orange. 



IV. RESULTS AND DISCUSSIONS 
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For C3:3 structure [Figure [TJa)], we find that the con- 
figuration with opposite spin orientation between adja- 
cent edge C atoms is energeticahy favored in comparison 
with the configuration with parahel spins between edge 
atoms of two sublattices represented by the blue and red 
color in Figure [H The total energy calculations reveal 
quite large magnetic interaction energies. For instance, 
the energy difference between ferromagnetic (FM) spin- 
polarized and paramagnetic (PM) state is found to be 
0.129 eV per edge atom. The spin configuration is fur- 
ther stabilized by 0.093 eV per edge atom as a result of 
the antiferromagnetic (AF) coupling between neighbor- 
ing atoms with magnetic moment of 0.48 /i^ per edge 
atom for each spin on each sublattice with opposite ori- 
entation. The magnetic moment is slightly larger than 
that of graphene nanoribbons which is ~0.43 fiB^^ • 

We now discuss the case of a 3-ring defect Ci2:i2 GNM 
[Fig. IHc)]. The total energy calculations show that the 
ground state is AF with a magnetic moment of 0.45 113 
per edge atom for each spin on each sublattice with op- 
posite orientation. The FM state is lower by 0.127 eV per 
edge atom compared to PM state, and the spin configura- 
tion is further stabilized by 0.107 eV per edge atom as a 
result of the AF coupling between neighboring atoms on 
different sublattice with opposite spin orientations [see 
TABLE [I]. 

The 1-ring and 3-ring defect C3:3 and Ci2:i2 structures 
considered above present no net permanent magnetic mo- 
ment, since the spin-polarized edge atoms appear in pair 
with opposite orientations, resulting in AF ground state 
with balanced spin- up and spin-down sublattices. How- 
ever, in a view of spintronic applications, it would be 
much more interesting to find the GNM structures with 



FIG. 2: H-passivated GNMs with triangular shapes: (a) Tq.^^ 
(b)Tfo:i2, (c)Tf5:i8, (d) Tfi^ss; with circled shape (e) C?2:i2] 
with rhombic shape (f) R24:24; with sector shape (g) Sf|:2i; 
and with pentagon shape (h) Pi8:2i- The corresponding net 
magnetic moments for each structure are also indicated. 



TABLE L The number of removed atoms on A and B sites and 
their difference A as- magnetic moment M(/is), total ener- 
gies E (in eV) for ferrimagnetic(FMi), antiferromagnetic (AF) 
and nonmagnetic (NM) states and defect formation energies 
E/ (in eV) for different GNM structures. The non-passivated 
C3:3, C6:6,Ci2:i2, and Tg:? and passivated Cf|:i2 GNMs are 
calculated with 6x6 unit cell, others are calculated with 8x8 
unit cell. 
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FIG. 3: (a) Total magnetic moment (/is/cell) (left) and spin- 
splitting (right) as a function of Aab for various GNM ge- 
ometries, where P'^-.iq is transformed from a pentagon struc- 
ture (Pi8:2i) by adding 2 A atoms to the 2 opened hexagons. 
The result of the Lieb's theorem prediction is also given for 
comparison. The even number of A+B structures Tf^.12 and 
T2i:25 are shown in red color to indicate the well agreement 
with Lieb's theorem prediction, (b) Energy difference between 
ferrimagnetic and paramagnetic states. 



nonzero net magnetic moment. This can be actually done 
by building unbalanced sublattice [Fig. [TJd)]. We found 
that the ground state of this unbalanced-defect triangular 
structure Tq.^j turns out to be ferrimagnetic (FMi) with 
total net moment of 3.987 /hb per unit cell. This moment 
originates from each edge atom's dangling bond (a-bond) 
with spin moment of 1 /i^ providing the magnetic mo- 
ment of 6(red)-3(blue)=3/is in addition to contribution 
from TT-bond equal to 1 fiB according to total number 
difference between atoms on A and B sublattice Aab- 

Dangling bonds at the edge C atoms of non-passivated 
GNM are strongly chemically reactive^, which leads to 
hole formation energy of pure GNM higher than 1 eV/C 
[see TABLE |T] . Thus edge C atoms are likely to be pas- 
sivated by light elements. We used hydrogen for passi- 
vation of edge C atoms and considered basic geometrical 
GNM shapes (Fig. [2j) for which the GNM hole formation 
energy is dramatically decreased due to passivation of 
dangling bonds [see TABLE |I]. For triangular holes, one 
can see that the formation energy decays as a function 
of hole dimension. At the same time, the magnetic mo- 
ment increases and is roughly proportional to the GNM 



hole size. When the triangular hole size is increased, 
one observes that the net moment gets bigger [Fig. [2jb,c 
and d)] and reaches 3.62 /i^ for the biggest hole shown 
in Fig. [2fd). It is interesting to note that Lieb's theorem 
was originally formulated for even A+B number of atoms, 
and indeed the obtained values for and T2^.25 fol- 

low Lieb's theorem predictions. However, there is an ex- 
ception for non-triangular case of (8^9.21) with total sum 
of A and B being even, which is not well accounted by 
Lieb's theorem. In addition, one can see from TABLE 
H] that the formation energy values of triangular GNMs 
decrease as a function of hole size and are comparable to 
those of the non-magnetic configurations. 

In Fig. [3l^a) we summarize aforementioned results in- 
cluding the calculated net magnetic moments for circu- 
lar (0^2:12)7 rhombic (R24:24)5 sector (829:21) and pen- 
tagon (Pi|:2i) GNM shapes represented in Fig. [2fe)-(h), 
respectively. In addition, the curves contain the net 
magnetic moment values for alternative pentagon shape 
GNM, 

Pif:i9 5 obtained from Pi|.2i by adding 2 A sites 
to complete 2 hexagons in upper left and upper right 6- 
rings in Fig. Efh). Even though the overall trend of the 
calculated values qualitatively follows the Lieb's theory, 
differences are observed, and first-principles calculations 
do not always correspond to Aab , even for the case with 
A+B is even of 8^9.21 as we have already mentioned. In 
fact, we can ascribe the structures with odd number of 
A+B atoms, i.e. Tgfy, Pi^i9, T^.^^g, Pf|.2i as well as sec- 
tor shape 8^9-21 GNMs to intermediate regime between 
nonmagnetic and highly magnetic regimes. This regime 
provides a root towards design of magnetic GNM super- 
meshes. It is worth to note that the possible mechanism 
for the deviation from the Lieb's theorem of the moment 
value for the sector shape GNM 8^9.21 compared to T 1^.12 
where A+B is even for both, could be attributed to larger 
amount of armchair edges (not favorable for moment for- 
mation) in 8^9.21 structure. 

To further elucidate the origin of magnetism in GNM 
structures, we compute the total and partial density of 
states (D08) for considered GNMs. Fig.|ll^a)-(d) give the 
total D08 for triangular shape GNMs shown in Fig.[2ja)- 
(d), respectively. The exchange splitting between major- 
ity and minority spins mainly originates from p^ orbit als, 
as clearly seen from Fig. IH^c') where the projected den- 
sity of states (pD08) on edge atoms for T^.^g is plotted. 
More interestingly, exchange splitting and energy differ- 
ences between FM and PM states also increase with A^^ 
following the same trend as the net magnetic moment [see 
Fig.[3l^a) and (b)], reaching values of 0.5 eV and 0.12 eV, 
respectively. These large exchange splitting values sug- 
gest that the magnetism could be preserved at room tem- 
perature which look very promising for room temperature 
graphene spintronics developments. Density of states of 
structures like T^.;l2 provides with the evidence for the 
presence of localized electrons on the zigzag edges similar 
to the perfect GNR^^. The smearing of the p^ peak and 
increase of the number of states at the Fermi level (met- 
allization) suggest the electron delocalization on edges 
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FIG. 4: Density of states for triangular GNMs of (a) T^fy, (b) 

Tfoa2, (C) Tfs.is, (d) (e) C|^3, (f) R?4:24, (g)P'lf 19, (h) 

Si9:2i5 and (i) Pi8:2i ^he peaks around Fermi level are marked 
with arrows. It can be seen that only state contributes to 
the moment from (c^) projected density of states of one edge 
atom in the Tf^^ig GNM. 



and switch from AFM to FM configuration similar to 
the case of partially oxidized graphite edges^^. We have 
plotted the spin density figures for the localized and de- 
localized cases for illustrate dramatic changes in localiza- 
tion of electrons with vanishing of pseudogap in density 
of states (see Fig. [SJ.The cause of such delocalization of 
the unpaired electrons on the edges of GNM holes can 
be attributed to combination of sublattice degeneracy 
breaking and deviation from the perfect shapes of the 
graphene nanoribbon. 
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FIG. 5: Spin density (/is/A^) distribution for the two types 
of graphene nanomesh with localized (upper panel) and delo- 
calized (lower panel) unpaired electrons. 



V. CONCLUSIONS 

In conclusion, electronic and magnetic properties have 
been explored in GNM with different geometries using 
first-principles calculations. For balanced non-passivated 
GNMs, the ground state was found to be either param- 
agnetic or antiferromagnetic. In situation of sublattice 
degeneracy breaking, ferrimagnetic ground states were 
obtained. The hydrogen-passivated GNMs were found 
to be strongly sensitive to the GNM size and shape, 
with magnetic moments deviation from the Lieb's the- 
orem trend caused by the delocalization of the unpaired 
electrons on the zigzag edges. The calculations of the for- 
mation energy provide with the evidence for the struc- 
tural stability and high probability of the formation of 
ferrimagnetic structures. Furthermore, three magnetic 
regimes are revealed: (i) highly magnetic GNM obey- 
ing Lieb's theorem ; (ii) quenched magnetic state due to 
complete chemical bond reconstruction; and (iii) interme- 
diate regime providing a possible root towards design of 
magnetic GNM supermeshes. These results demonstrate 
that a turn from zero dimensional graphene nanoflakes 
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throughout one- dimensional GNR with zigzag edges to 
GNM breaks locahzation of unpaired electrons at zigzag 
edges and provides a deviation from the Lieb's theorem 
trend. Such delocalization of the electrons allows switch- 
ing of the ground state of system from antiferromagnetic 
narrow gap insulator discussed for GNRs to ferromag- 
netic or nonmagnetic metal. These results combined with 
obtained large values of the exchange splitting (increasing 
with A^^) pinpoint promising perspectives for develop- 
ing room-temperature graphene spintronics. 
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